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I  INTRODUCTION 

This  is  the  seventh  quarterly  technical  sunsaary  report  issued 
by  SDL  and  covers  the  period  October  through  December  1966,  Work 
previously  completed  or  work  currently  in  progress  is  mentioned 
only  ss  it  relates  to  analyses  coexisted  during  this  reporting 
period. 

Analyses  performed,  for  which  results  have  been  reported 
are  discussed  in  Section  1 1  under  descriptive  headings,  section 
III  contains  a  discussion  of  the  support  and  service  tasks  per¬ 
formed  for  in-house  projects  and  for  other  VELA  UNIFORM  participants. 

Appendix  A  is  a  listing  of  those  organisations  receiving  SLL 
data  services  during  this  period;  and  Appendix  B  contains  selected 
SDL  reports  which  are  representative  of  the  types  of  analyses  made 
by  the  Seismic  Data  Laboratory. 

II.  WORK  COMPLETED 

A*  Finite  Fourier  Transform  Theory  and  Itw  Application  to  the 
Computation  of  Convolutions,  Correlations,  and  Spectra 

Practical  and  computational  aspects  of  the  theory  of  Fourier 
Transforms  have  been  examined  in  connection  with  various  SDL  research 
analyses.  These  efforts  have  resulted  in  a  set  of  programs  for  per¬ 
forming  operations  on  time  series  based  on  the  Cooley-Tukey  (Refer¬ 
ences  1,2)  hyper-rapid  Fourier  transform  method.  Using  this  method, 
computations  on  seismic  array  data  such  as  the  calculation  of  con¬ 
volutions,  correlations,  spectra,  and  digital  filters  have  been 
speeded  up  by  factors  of  three  or  four  and  sometimes  even  ten.  Fol¬ 
lowing  is  a  brief  description  of  the  procedures  followed  and  the 
results  obtained  from  these  procedures. 
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1,  Finite  end  Discrete  Fourier  Transform# 

In  the  cete  of  continuous  date  of  infinite  length,  the 
Fourier  transform  pair  ia  the  direct  transform  and  the  inverae  trans- 
forai,  Sometimes  the  direct  transform  ia  written  with  a  factor  of  1 
in  front  on  the  integral  and  the  inverae  with  a  factor  of  1/2-  . 
Ouantitiea  of  intereat,  auch  aa  apectra,  etc.,  involve  magnitudes  or 
squares  of  one  transform  and  the  factor  suit  be  inserted  or  taken 
out,  dapending  on  which  definition  is  used,  to  preserve  the  proper 
dimension#. 

TVo  drawbacks  of  these  definitions  for  digital  confuta¬ 
tion#  are  apparent:  First,  the  integrals  must  be  approxiausted  by 
sums  in  the  digital  computer,  v*jich  implies  that  both  transforms  in¬ 
volve  sampled  variables.  Second,  the  infinite  limits  on  the  sums 
art  impossible  in  practice.  Clearly,  these  sums  must  be  truncated, 
as  they  do  not  in  general  converge  over  a  finite  interval.  As  a 
result,  Fourier  transforms  as  such  are  never  really  computed  by  a 
digital  computer  Instead,  the  complex  samples  of  a  direct  trans¬ 
form  are  approximated  by  the  cosine  and  sine  coefficients  of  Fourier 
series  representation  of  the  input  data. 

If  N  samples  of  the  data  are  taken  at  equally  spaced 
intervals  6t  -  T/R,  the  integrals  become  sums,  and  the  frequency 
sum  goes  from  DC  to  the  folding  frequency.  it  was  found  that  a 
great  deal  of  symmetry  between  the  two  transforms  could  be  pre¬ 
served  if  the  sum  is  summed  up  to  M-l,  Redundant  points  in  the 
spectrum  are  included  (since  the  transforms  are  periodic)  but  the 
computational  procedures  are  simplified. 

w*a  shown  t  •  the  set  of  direct  Fourier  trans¬ 
form  points,  between  DC  and  the  folding  frequency,  contained  the 
same  amount  of  information  as  the  real  data  series,  which  suggested 
that  tne  existence  of  one  transform  should  imply  the  existence  of 
the  other.  » 


2.  Two-and-Threc-Diraensional  Fourier  Transforms 

Two-  and  three-dimensiona 1  direct  Fourier  transforms 
are  seen  to  be 
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By  separating  and  breaking  up  these  equations,  it  was 
calculated  that  +  iS^  one-dimensional  transforms  are  required  to 
conpute  the  single  two-dimensional  transform,  and  that  one¬ 

dimensional  transforms  and  two-dimensional  transforms  are  needed 
to  compute  the  single  three-dimensional  transform. 

3.  Speed-Up  of  Transform  Computing  Time 

It  was  shown  that  the  process  of  transforming  is 
equivalent  to  matrix  multiplication  by  a  matrix  W,  which  preserves 
"length"  between  two  domains.  The  Cooley-Tukey  method  factors  the 
w  matrix,  it  its  order  is  a  power  of  two,  into  L  ♦  1  sparse  matrices 
where  L  is  th*  power  of  two.  Multiplying  L  ♦  1  times  by  these  apart 
matrices  can  in  some  cases  reduce  the  computing  time  by  many  tens  of 
times , 


4.  High-Speed  Cor relations  and  Convolutions 

By  computing  Fourier  transforms  with  a  finite  Fourier 
series-like  method  an  important  condition  -*  put  on  the  time  series. 
As  in  regular  Fourier  series  the  input  is  assumed  to  be  periodic 
with  period  T  and  the  integrals  or  sums  are  computed  over  a  single 
period.  There  iu  also  the  effect  of  cutting  off  the  spectrum  at 
the  folding  frequency.  Sines  and  cosines  of  finite  wavelength  will 
repeat  again  outside  the  region  of  interest.  This  fact  in  itself  is 
not  bothersome  but  becomes  a  serious  complication  in  the  computation 
of  convolutions  and  correlations.  Convolutions  and  correlations  as 
usually  computed  assume  the  time  series  to  be  zero  outside  the  region 
of  interest.  Therefore,  the  integrals  or  sums  in  computing  them  are 
sunsned  out  only  over  the  non-zero  terms.  When  multiplying  together 
two  finite  Fourier  transforms  (or  the  complex  conjugate  of  one  times 
the  other)  the  periodicity  of  the  time  series  means  that  elements 
whicn  have  been  shifted  past  the  end  of  a  period  reappear  at  the 
beginning.  This  process  is  called  circular  convolution  or  correla¬ 
tion  and  its  effects  are  unavoidable  when  straightforwardly  comput¬ 
ing  lagged  products  with  finite  Fourier  transforms. 

Circular  convolution  is  written: 

T-l 

RC,(t)  -  E  x  (t)  x  (t+f  )  (3) 

*■3  -r»0  1  J 

where  x  (t  +  T)  -  x  (t)  for  all  m 
m  m 

It  was  shown  that  this  kind  of  correlation  is  equal 
to  the  transform  of  the  absolute  product  of  the  two  finite  trans¬ 
forms. 
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On  the  other  hand  the  transient  correlation  for 
positive  lags  is  defined  by  the  following: 

T  T-l-t 

R  (t)  =  r  x  (t)  X  (t  +  T)  (4) 

where  the  upper  limit  on  the  sum  simulates  the  desired  zeros  in 

the  time  series  outside  the  region  of  interest.  The  finite 
,  t 

Fourier  transform  of  this  R  is  thus  not  the  product  of  the  two 
individual  transforms  However,  by  filling  zeros  into  the  second 
half  of  each  data  series  and  computing  their  transforms  out  to 
twice  their  actual  length,  a  good  estimate  of  the  spectrum  may  be 
obtained.  In  addition,  the  negative  lags  in  the  correlation  appear, 
thus  giving  a  more  mathematically  satisfying  result. 

Transient  correlations  for  100%  lags  were  shown  to  be 
computed  by  forming  the  absolute  product  of  two  transforms,  each 
computed  out  to  twice  the  length  of  the  original  data  series  with 
zeros  filled  into  the  second  halves 

Non-circular  or  transient  convolutions  were  also  com¬ 
puted  in  much  the  same  way,  except  that  the  transforms  had  to  be 
computed  out  to  a  length  equal  to  the  sum  of  the  lengths  of  the 
time  series  and  the  filter,  with  the  appropriate  number  of  zeros 
filled  into  each.  The  convolution  theorem  was  proved  in  the  same 
fashion  and  were  computed  by  forming  the  product  of  the  two  trans¬ 
forms,  each  computed  out  to  a  length  equal  to  their  sum  with  zeros 
filled  into  the  extra  lengths. 
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B‘  Mylgjah  Wave  Rejection  by  Ontinv  m  Filterina  of  Vernal 
Arrays  — ~ - — 

Vertical  arrays  offer  some  intriguing  advantages  over  other 
arrangements  of  seismometers  Assuming  the  noise  to  be  composed  of 
Rayleigh  waves  and  perhaps  some  mantle-propagating  P-waves,  and 
assuming  the  signals  to  be  teleseismic  P-waves,  the  noise  and  signals 

will  be  recorded  by  a  vertical  array  with  some  rather  special  charac¬ 
teristics. 


We  can  imagine  that  each  mode  of  Rayleigh  wave  noise  pos¬ 
sesses  a  random  character  common  to  noise  functions  in  general. 
However,  for  each  Rayleigh  wave  '.ode  the  response  versus  depth  and 
frequency  relative  to  its  response  on  the  surface  is  predictable. 
Because  the  relative  depth  variation  is  predictable,  the  vertical 
array  can  be  summed  to  cancel  the  Rayleigh  modes,  m  contrast, 
the  noise  over  a  surface  array  is  more  unpredictable  and  more 
varii ble  and,  therefore,  more  difficult  to  cancel  by  array  summation. 

The  purpose  of  this  analysis  was  to  test  the  theory  of 
applying  maximum  likelihood  filters  in  vertical  arrays  in  order  to 
provide  undistorted  estimates  of  the  signal  or  the  various  Rayleigh 
modes  with  the  other  modes  cancelled  out. 

A  vertical  array  of  five  vertical  component  seismometers 
which  existed  in  a  well  at  the  Uinta  Basin  Seismological  Observatory 
(UBSO)  was  used  for  testing  this  theory.  The  Rayleigh  dispersion 
curves  for  this  well  indicated  that  the  fundamental,  first,  and 
second  higher  Rayleigh  modes  would  all  be  supported  over  the  signal 
frequency  range  of  0.5  to  2.0  cps  Since  we  wanted  more  recording 
levels  than  modes,  these  three  Rayleigh  modes  plus  a  signal  mode 

were  the  only  ones  considered  to  exist  in  our  synthetic  modeling 
of  this  well* 
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The  optimum  filter,  G,  solutions  are  the  least  squares 
inverse  to  the  Rayleigh  signal  filter  matrix,  H.  Thus,  the  mitrix 
products  of  G  and  H  gives  the  identity  matrix.  In  our  test  model 
we  assumed  the  signal  existed  with  the  same  size  and  amplitude  on 
all  traces  We  derived  optimum  filter  solutions  only  over  the 
signal  range  of  0.5  to  2.0  cps.  For  the  frequencies  below  0.5, 
we  smoothed  the  frequency  responses  to  zero  according  to  a  sinusoidal 
1J *  f unction  for  both  the  H  and  G  frequency  responses.  We  smoothed 
these  responses  to  zero  in  a  similar  way  over  the  range  from  2  to 
5  cps.  Also,  the  orthogonality  of  our  optimum  filter  solutions  were 
maintained  over  the  entire  frequency  range. 

In  order  to  demonstrate  that  our  optimum  filter  solutions 
were  correct,  we  created  some  synthetic  data  which  conformed  to  the 
well  locr  analysis  we  assumed  for  the  UBO  well.  These  Rayleigh  mode 
traces  were  combined  with  an  artificial  signal  to  produce  a  mixture 
of  signal  and  Rayleigh  modes  expected  from  this  well. 

Figure  1  shows  the  data  from  the  five  levels  of  the 
vertical  array  if  the  signal  plus  the  fundamental  and  the  first  and 
second  higher  Rayleigh  modes  have  been  added  into  the  vertical  array 
data.  The  outputs  for  the  signal  and  the  three  Rayleigh  modes  are 
reproduced  in  size  and  waveform  within  a  few  percent  of  the  expected 
results.  The  signal  trace  shows  better  than  a  20  db  improvement 
over  the  best  signal/noise  ratio  available  in  the  vertical  array. 

Thus,  the  optimum  filters  perform  as  expected  when  the  data  proper¬ 
ties  match  the  dispersion  analysis  of  the  well  log. 

These  solutions  were  then  applied  to  actual  noise  data 
recorded  at  the  vertical  array.  All  of  the  optimum  filters  used 
required  inputs  from  all  five  levels  of  noise  data.  However,  we 
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could  ask  for  estimates  of  only  a  single  mode  or  for  estimates  of 
several  modes  in  i.he  filter  outputs.  Th^  more  modes  a  set  of 
filters  estimate,  the  larger  the  filter  gains  will  be.  It  follows 
that  any  errors  resulting  from  incorrect  dispersion  analysis,  in¬ 
strument  calibrations,  or  non-parallel  layering  will  be  amplified 
much  more  in  these  filters  with  fewer  degrees  of  freedom.  Results 
showed  that  the  signal  estimate  outputs  gave  progressively  increasing 
low  frequency  errors  on  the  estimates  of  all  modes  as  the  degrees  of 
freedom  decreased  (i.e.,  as  the  number  of  outputs  increased). 

It  was  concluded  that  our  optimum  filter  solutions  have 
the  following  features: 

1.  These  filters  are  indeed  orthogonal  and  separate  the 
noi 3e  and  signal  modes  as  planned. 

2.  These  optimum  filters  extend  for  400  points  in  time 
since  the  frequency  interval  chosen  was  0.5  cpa. 

3.  The  gain  on  the  synthetic  examples  was  greater  than 
20db  with  solutions  restricted  to  the  0.5  -  2.0  frequency  range. 

More  resolution  in  frequency  could  increase  this  figure. 

4.  The  optimum  filters  are  zero  phase  shift  filters. 
Therefore,  they  are  algebraically  additive. 

5.  The  optimum  filters  become  larger  in  gains  (both 
positive  and  negative)  as  degrees  of  freedom  decrease. 

6.  The  optimum  filters  become  larger  in  gains  (both 
positive  anc  negative)  as  the  aperture  of  array  goes  down. 

Thus,  for  low  frequencies  the  solutions  tend  to  become  unstable. 

7.  Extra  modes  not  considered  and  errors  in  assumptions 
cause  errors  in  the  output.  Errors  in  assumptions  can  include  an 
incorrect  well  log,  non-parallel  layering  in  surrounding  medium, 
and  incorrect  calibrations  of  seismometers. 


8.  Extra  degrees  of  freedom  are  needed  to  cut  gains  of 
optimum  filters  and  make  optimum  solutions  more  tolerant  of  errors. 

9,  Extra  degrees  of  freedom  are  best  obtained  by  increas¬ 
ing  the  number  of  seismometers  in  the  vertical  array. 

10.  More  stable  solutions  (i.e.,  optimum  filters  with 
lower  gain)  will  be  obtained  rrom  the  deeper  vertical  arrays  which 
have  the  seismometers  distributed  rather  uniformly  throughout  the 
array. 

C.  LASA  Travel  -Time  Data  at  the  SDL 

In  connection  with  ite  many  LASA  analyses,  the  SDL  is 
accumulating  relative  travel-times  at  the  Montana  Large  Aperture 
Seismic  Array  (LASA)  with  a  view  towards  computing  travel-time 
anomalies  at  the  21  subarrays.  In  the  belief  that  these  raw  data 
are  in  demand  and  may  be  of  use  to  the  seismic  community  in  general, 
a  report  was  issued  which  contains  LASA  travel-time  data  for  approx¬ 
imately  400  events,  as  read  from  LASA  films. 

Figure  2  shows  the  type  of  data  which  has  been  computed 
and  the  following  numbered  explanation,  corresponding  to  the  circled 
number  on  the  figure,  describes  the  method  of  presentation] 

1.  Arbitrary  region  name 

2.  Direction  of  approach 

3.  Distance  range  for  events  included  in  region 

4.  Azimuth  range  for  events  included  in  region 

5.  Event  date  and  name 

6.  PDE  latitude 

7.  PDE  longitude 

8.  PDE  depth 


9.  PDE  origin  time 
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Figure  2.  Type  of  LASA  Travel-Time  Data  Computed  by  the  SDL 
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10.  Arrival  time  in  seconds  (an  arrival  time  of  "0" 
indicated  no  reading  made  at  that  subarray. 


First  row,  left  to  rights 

Bl, 

B2, 

B3, 

B4, 

Cl, 

C2, 

C3 

Second  row  : 

C4, 

Dl, 

D2, 

D3, 

D4, 

El, 

E2 

Third  row  : 

E3, 

E4, 

FI, 

F2 , 

F3, 

F4, 

AO 

11.  Arrival  time,  hour  and  minute  at  LASA 
D.  Energy  Fluctuations  in  Seismic  Noise 

In  recent  years,  statistical  methods  have  been  successfully 
used  to  predict  the  change  in  ambient  seismic  noise  power  sensed  by 
burying  a  seismometer  beneath  the  surface  of  the  earth.  Observa¬ 
tions  of  seismic  noise  has  indicated  a  large  standing  wave  or  iso¬ 
tropic  component  composed  of  the  admixture  of  many  propagation  modes. 
Statistical  theories  such  as  eguipartition  of  energy  have  been  reason¬ 
ably  successful  for  deriving  the  excitation  of  the  propagation  modes. 

Since  such  models  are  classically  described  by  waves  from  a 
zero-mean  Gaussian  population,  the  purpose  of  this  study  was  to 
tentatively  assume  that  the  seismic  noise  is  Gaussian  and  to  try  to 
reject  the  hypothesis  by  measuring  the  relative  frequency  of  occur¬ 
rence  of  specific  particle  energies  seguentiallv  observed  on  samples 
of  seismic  noise.  As  control  for  the  study,  the  same  frequency  of 
occurrence  observations  were  made  of  a  thermal  noise  sample  recorded 
on  tape  from  a  Gaussian  noise  generator.  The  energy  envelope  of 
Gaussian  noise  is  theoretically  described  by  a  Boltzman  or  exponential 
probability  function,  thus  the  fit  of  a  straight  line  to  the  log  of 
the  relative  frequency  of  occurrence  is  a  test  of  the  Gaussian  hypoth¬ 
esis.  Two  tests  were  used;  one  was  the  chi-square  test  of  the  devia¬ 
tion  from  the  least  squares  straight  line;  the  other,  more  powerful 


test,  compared  the  variance  of  the  seismic  noise  and  that  of  the 
Gaussian  noise  generator. 

Three  hours  of  seismic  noise  was  selected  which  appeared 
to  be  typical  of  the  normal  ambient  background.  The  noise  was 
measured  at  a  depth  of  7,452  feet  in  a  deepwell  near  Apache, 

Oklahoma  The  near  surface  layering  i«  described  as  high  velocity 
limestones  of  approximately  6  km/aec  overlying  an  igneous  basement 
complex  of  velocities  in  the  neighborhood  of  5  kra/sec. 

As  shown  on  Figure  3  for  seismic  noise,  the  relative  fre¬ 
quency  of  occurrence  is  given  by  the  vertical  axis,  and  the  class 
interval  in  order  of  increasing  energy  is  shown  on  the  norixontal 
axis,  on  each  plot,  the  center  frequency  and  effective  bandwidth 
of  the  filter  is  labeled.  On  Figure  4,  the  same  information  is 
shown  for  the  Gaussian  noise  generator.  Comparing  Figure  3  with 
Figure  4,  the  scatter  from  a  linear  trend  is  the  same  overall  for 
both  sets  of  data,  seismic  and  Gaussian  generator. 

Based  on  the  chi-square  test,  it  is  more  than  99%  probable 
that  samples  at  all  of  the  frequencies  are  from  a  zero  mean  Gaussian 
population.  The  hypothesis  that  the  class  deviation  of  the  seismic 
noise  sample  and  the  thermal  noise  sample  are  from  the  same  random 
population  is  rejected  only  for  the  sample  at  1.4  cps,  suggesting 
that  a  non-Gaussian  component  may  be  included  with  the  noise  in 
this  band.  if  we  accept  the  Gaussian  hypothesis,  theoretical  in¬ 
terest  in  the  power  spectral  density  of  the  noise  may  bear  ultimately 
on  dissipation  mechanisms  and  may  possibly  lead  to  analysis  of 
geological  structure,  with  the  real  changes  in  the  power  spectrum 
m  a  region  depending  primarily  on  changes  in  structure.  Of  more 
pra  t-ical  interest  is  that  modern  literature  on  detection  and  fil¬ 
tering  is  roost  meaningful  in  the  context  of  seismic  signals  added 
to  Gaussian  noise. 


11 


E.  Shot  and  Earthquake  Analyses 

Standard  atatiatical  summary  reports  were  issued  on  the 
Nevada  Teat  Site  (NTS)  events  TAN  and  HALF  BJAK.  TAN,  detonated 
on  3  June  1966  in  a  tuff  geologic  medium,  had  an  average  magnitude 
of  5.56  ♦  0.49;  ahort-period  aignala  were  recorded  by  22  stations 
and  long-period  aignala  were  recorded  by  21  stations.  HAL.  BEAK, 
detonated  on  30  June  1966  in  a  rhyolite  geologic  medium,  had  an 
average  magnitude  of  6.02  +  0.60,  with  ahort-period  and  long-period 
aignala  being  recorded  by  21  and  20  station*,  respectively. 

Figure  5,  showing  the  Lg  amplitudes  for  TAN  and  HALF  EEAK,  are 
illustrative  of  the  atatiatica  summarized  from  these  analysea. 

Ocher  projects  completed  included  the  following: 

-  signal  albums  (photographic  enlargements  of  t  16  mm 
and  35  mm  film  records  for  the  LRSM  and  observatory 
stations)  from  5  events. 

-  analysis  of  575  events  for  a  special  long-period  study 

-  film  analysis  of  CPSO  noise 

-  special  LASA  PcP  study 

F  vertical  Array  Teleseismic  Signal  Measurements 

The  concept  of  a  vertical  array  is  to  record  simultaneously 
on  several  transducers  stacked  in  a  deepwell  3  km  or  more  in  depth. 
The  purpose  of  this  study  wss  to  investigate  the  possibility  of 
reducing  near  surface  reverberations  due  to  geological  effects  near 
the  vertical  array  receivers,  based  on  procedures  established  over 
a  pei iod  of  time  in  reflection  seismology,  where  vertical  stacking 
of  an  explosion  source  is  used  to  remove  distortion  of  near  surface 
layering,  mainly  the  ghost  reflection  from  the  surface.  By  reci¬ 
procity,  the  same  algorithms  developed  for  stacked  sources  can  be 
applied  to  stacked  receivers  by  reversing  the  sign  of  the  reflection 
coefficients. 
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Figure  5.  Lg  Amplitudes  for  Two  NTS  Nuclear  Explosions 


The  data  used  were  three  events  recorded  it  a  deepwell, 
each  having  different  focal  depths  and  different  apparent  fre¬ 
quencies. 

F«=fore  developing  the  deghosting  technique,  the  deepwell 
trace  was  simulated  from  a  surface  seismogram,  if  a  deepwell  trace 
can  be  simulated  by  adding  only  an  echo  from  the  surface,  it  seems 
entirely  plausible  that  by  removing  the  ohost  reflection  from  the 

deepwell  trace,  the  up-going  P-pulses  can  be  made  to  appear  similar 
at  each  depth. 

Two  deghosting  techniques  were  developed  as  follows: 

1*  When  a  surface  or  near  surface  trace  is  available,  the 
deghosted  trace  can  be  constructed  by  simply  shifting  the  surface 
trace  by  K,  multiplying  a  Jja,  and  subtracting  this  resultant  from 
either  the  deepwell  trace  (phased  vertical  array)  or  the  deepwell 
trace  shifted  by  K/2  (non-phased  vertical  array) . 

2.  It  would  seem  more  plausible  to  develop  a  deghosting 
process  that  removes  the  echo  without  the  use  of  a  surface  seismo¬ 
gram,  a  practical  technique  that  would  figuratively  push  the  ghost 
off  the  end  of  the  deepwell  trace. 

Given  the  reflection  coefficient  a  ,  and  the  deepwell 

tia  :e  Yi  x  the  phased  deghosted  trace  is  formed  as  follows:  First  a 
new  trace  Z,  .  is  defined  to  be 

Zi>  1  *  Yi/1  ~  X(i  -  NT^  *  ai, 

where  is  equal  to  the  product  of  the  sampling  rate  and  the  echo 
time  delay  fK) .  This  above  equation  simply  shifts  and  invv rts  the 
ghost  to  a  point  farther  down  the  seismic  record  than  it  originally 

was.  This  iterative  process  is  continued  until  the  ghost  reflection 
is  pushed  off  the  seismic  record. 


Results  tended  to  show  that  vertical  array  processing, 
based  on  previous  experience  with  stacked  sources,  appears  to  br 
an  effective  means  of  reducing  the  reverberations  and  complexity 
of  the  coda  caused  by  near  surface  reverberations  at  the  receiver. 

In  Figure  6,  the  first  two  traces  show  the  raw  data  measurements 
oud  its  correlation  traces,  the  next  two  traces  show  the  phased 
deepwell  measurements  deghosted  by  subtraction  of  the  surface 
measurement,  and  the  last  two  traces  show  the  phased  deepwell 
measurements  deghosted  by  means  of  an  inverse  operator  which  removes 
the  surface  echo.  Both  deghosting  methods  performed  adequately. 

The  correlation  trace  for  the  inverse  operation  had  a  slightly  cleaner 
coda  and  this  method  has  the  advantage  of  not  requiring  a  quiet  sur¬ 
face  trace  and  exact  control  of  instrument  gain.  The  correlation 
traces  used  for  detecting  the  upgoing  P-pulses  impose  the  strongest 
possible  requirement  that  che  signal  be  fixed  jointly  on  all  channels. 
III.  SUPPORT  AND  SERVICE  TASKS 

A.  VELA  UNIFORM  Data  Services 

As  part  of  the  contract  work-statement,  the  SDL  provided 
one  or  more  of  the  following  support  and  service  functions  for  VSC 
and  other  VELA  participants: 

-  copies  of  16  and  35  mm  film 

playouts  of  earthquakes  and  special  events 

-  copies  of  existing  composite  analog  tapes 
composite  analog  tapes  of  special  events 

-  use  of  1604  computer  for  checking  out  new  programs 
or  running  production  programs 

-  copies  of  digital  programs 

digitized  data  in  standard  formats  or  special  formats 
for  use  on  computers  other  than  the  1604 
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Figure  6.  Deghosting  Techniques  for  Reducing 
Near  Surface  Reverberations 


-  running  SDL  production  programs,  such  as  power 
spectral  density  and  ar.ay  processing,  on 
specified  data 

-  digital  x-y  plots  of  power  spectra  or  digitized 
data 

-  signal  reproduction  booklets 

-  copies  of  MIT  Geophysics  Program  Set  II 

-  space  for  visiting  scientists  utilizing  SDL  fa¬ 
cilities  to  study  data  and  exchange  information 
with  SDL  personnel. 

During  this  report  period,  80  such  projects,  including  those 
for  the  VELA  Seismological  Center,  were  completed  and  the  18  organ¬ 
izations  receiving  these  services  are  listed  in  Appendix  B. 

B.  Data  Library 

The  Data  Library  contains  approximately  7,484  digitized 
seismograms,  165  digital  computer  programs  and  296  composite  analog 
magnetic  tapes,  all  available  for  use  by  the  VELA  UNIFORM  program. 

The  following  additions  were  made  during  this  report 

period: 

1.  Digital  Seismograms  -  750  including 

-  8  earthquakes  recorded  at  one  station 

-  2  earthquakes  recorded  at  24  stations 

-  6  earthquakes  recorded  at  85  stations 

-  long  period  data  from  7  earthquakes 

2.  LASA  Data  -  180  digital  tapes  and  210  16mm  films 

-  there  are  a  total  of  703  LASA  digital  tapes  and 
354  16mm  films  in  the  library.  There  is  also 
a  master  calibration  tape  which  containe  the 
magnification  (digital  counts  per  millicron)  of 
each  sensor  for  every  subarray.  These  magnifi¬ 
cations  have  been  computed  for  all  calibration 
tapes  currently  in-house.  As  each  new  calibration 
is  received,  it  is  routinely  run  through  the 
program  CALIBR  and  added  to  the  master  tape. 


3.  Digital  Programs  -  15  including: 

SVARRAY  -  To  process  vertical  array  data  aligned  for 

first  P-motion. 

FASTFIL  -  To  broadband  pass  filter  seismograms. 

GOSTSEIS-  To  process  vertical  array  data.  This  package 
removes  ghost  reflection,  forms  two  sum  traces  -  one  the  deepwell, 
the  other  the  deghosted  deepwell  and  computes  a  recursive  corre¬ 
lation  trace  using  the  deghosted  deepwell  data.  Plots  are  obtained 
of  each  deepwell  and  deghosted  deepwell  seismogram,  the  surface 
trace,  the  two  above  sum  traces  and  the  correlation  trace. 

DEGSEIS  -  To  deghost  vertical  array  seismograms  given 
an  echo-time  and  reflection  coefficient  for  each  seismogram.  In 
addition,  a  sum  trace  of  the  deepwell  traces,  and  a  sum  trace  of 
the  aligned-deghosted-deepwell  traces  are  computed.  Also,  the 
component  which  is  jointly  correlated  on  all  aligned-deghosted- 
deepwell  traces  is  computed.  Furthermore,  this  program  plots  all 
deepwell  and  deghosted  deepwell  traces,  the  pair  of  sum  traces, 
and  the  correlation  trace. 

PGOSTSEI  -  Ibid  "Vertical  Array  Processor  Package  I  but 

uses  subroutine  SDEGGSTE  instead  of  DEGOSTE. 

SDEGSEI  -  To  remove  ghost  reflection  on  vertical  array 
seismograms  aligned  for  first  P-motion,  given  only  the  echo-times 
and  reflection  coefficients.  In  addition,  sum  traces  of  both  the 
aligned  deepwell  and  aligned-deghosted-deepwell  traces  are  computed. 
Also,  the  component  which  is  jointly  correlated  on  all  deghosted 
reismograms  is  computed.  Furthermore,  plots  are  obtained  of  all 
deepwell  and  deghosted  deepwell  traces,  the  surface  trace,  the  two 
above  sum  traces,  and  the  correlation  trace. 

VARRAY  -  To  process  vertical  array  data. 

RESPONSE  -  Tc  calculate  and  plot  in  k-space  the  response 
of  an  array  of  seismometers  with  respect  to  the  response  of  the 
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same  array  to  waves  wi.th  infinite  phase  velocity  (i.e.,  number  of 
db  down  from  zero) . 

ABML  -  To  find  the  maximum  absolute  value  of  o  fixed  or 
floating  point  number  in  an  array  and  the  FORTRAN  subscript  of 
the  arras  element . 

ROOTMNSQ  -  Using  standard  library  seismogram  tapes  and 
LASA  calibrition  tapes,  this  program  will  compute  the  mean  over 
any  portion  >f  any  cl  annel  for  a  particular  seismogram,  demagnify, 
filter  or  not  filter  (SUL  and  L/L  coefficients  available),  and  then 
compute  the  IMS  value  over  any  length  of  the  channel  desired. 

RICFIL3  -  To  perform  one-pass  shift  recursive  filtering 
of  seismic  cata  given  the  Q-value,  center  frequency,  and  sampling 
rate . 

REVERSE  -  To  time  reverse  or  flip  seismic  data. 

COLYTUKY  -  To  compute  Auto-spectra,  Cross-spectra,  Auto¬ 
correlations  and  Cross-correlations  by  using  the  COO'EY-TUKEY 
method . 

CNTUR4  -  Writes  printer  plot  of  the  contours  of  a  two- 
dimensional  array.  The  contours  are  formed  by  shading  the  area 
between  level  21  and  level  21+1.  This  display  is  superimposed 
upon  either  a  grid  work  or  a  border  with  tick -marks.  The  contour 
levels  may  be  specified  individually,  or  by  the  increment  between 
levels.  In  addition,  provision  has  been  made  to  allow  the  ori¬ 
entation  of  the  array  to  be  altered  by  either  exchanging  the  rows 
and  columns  or  reversing  either  the  rows  or  the  columns. 

GRTZEL  -  To  compute  the  Fourier  series  expansion  of  a 
rosi-valued  even  function,  or  the  real-valued  even  function  of  a 
real-valued  Fourier  series  expansion. 

4.  Analog  composite  Tapes  -  13  including: 

a.  Made  by  SDL 

-  9  special  composites 
-  17  - 


b.  Hade  by  Gaotech 

-  DUMONT 

-  HALF  BEAK 

-  PILE  DRIVER 

-  TAN 

C.  Data  Cotnpreaaion 

This  ia  a  continuing  routine  operation,  and  production 
la  maintained  at  the  level  needed  to  meet  the  requirements  of  the 
field  operation  (LRSH  and  U.S.  Observatories)  and  the  Seismic  Data 
Laboratory.  For  this  period,  3,767  tapes  were  compressed. 

D.  Equipment  Modifications 

The  A/D  system  has  been  modified  so  that  the  digitizing 
rate  can  be  synchronized  to  VELA  Time  Code,  thus  minimizing  speed 
differences  between  analog  recordings  from  different  stations. 

The  Pace  231-R  analog  computer  has  been  modified  to 
provide  additional  gain  input  to  each  amplifier  so  that  both  do 
not  have  to  be  utilized  when  increasing  gain. 

E.  Automated  Bulletin  Process 

July  and  August  1966  LRSM  and  Observatory  bulletins  vere 
processed  duving  this  report  period  and  forwarded  to  the  Geotech 
Division,  Teledyne  Industries,  Inc.,  for  checking  and  publication. 

F.  Special  Data  Processing 

A  1604  program  has  been  written  to  put  the  LASA  slow  mode 
data  tapes  into  the  SDL  library  format.  Either  the  Bhort  period 
center  instrument  data  or  the  long  period  three  component  data  can 
bt  processed. 

A  160-A  program  has  been  written  for  D  to  A  conversion 
of  the  LASA  slow  mode  tapes.  Thus  DAC  playouts  can  be  made 
directly  from  the  800  bpi  LASA  field  tapes. 
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ORGANIZATIONS  RECEIVINu_gPL  DATA  SERVICES 
October  1966  -  December  1966 

Brown  Engineering,  Incorporated 

Canadian  Pacific  Oil  and  Gaa  Company 
Continental  Oil  Company 

Earth  Sciences  Division,  Teledyne,  Inc. 

General  Atronics  Corporation 
Geotech  Division,  Teledyne,  Inc. 

Graduate  Research  Center  of  the  Southwest 

Massachusetts  Institute  of  Technology 
Mitre  Corporation 

Oregon  State  University 

Penn  State  University 
Princeton  University 

Rome  Air  Force  Base 

St.  Louis  University 

Texas  Instruments,  Incorporated 

University  of  Tasmania 

U.  S.  Coast  and  Geodetic  Survey 


Vitro  Corporation 
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Bradford,  J. 
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Chiburis,  E. 

Chiburis,  E. 

Chiburis,  E. 
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Measurement 


SEM£TEP-SPL  REPORTS 


C.  B. ,  Flinn,  E.  A.,*  Detection,  Analysis  and 
Interpretation  of  the  Telcseismic  Signal 
f r the  SALMON  Event;  2  April  1965 

C.  B. ,  Flinn,  E.  A.;  Perturbation  Theory  for 
the  Inversion  of  Body  Wave  Travel-Time  Data; 

5  April  1966 

.;  Estimation  of  Network  Capability;  20  January  1964 

. ;  Analysis  of  Variance  as  a  Method  for  Seismic 
Signal  Detection;  25  February  1965 

C. ;  Evaluation  of  the  TFSO  31  Element  Array; 

13  December  1963 

C. ,  Enochnon,  L.  D.  S  and  Thrall,  G,  P. *;  Seismic 
Partial  Coherecny;  28  April  1965 

C. ;  Enochson,  L.  D. Numerical  Experiments  with 
Partial  and  Multiple  Coherency;  29  October  1965 

F. ;  Ahner,  R.  O. ;  LASA  Travel-Time  Date  at  the 
SDL;  5  December  1966 

F.;  Dean,  W.  C. ;  Teleseismic  Signal  Alignment  at 
the  Tonto  Forest  Extended  Array;  15  October  1965 

F.  ,*  Relative  Travel-Time  Anomalies  at  LASA  and  the 
Location  of  Epicenters  Using  "SHIFT";  7  June  1966 

F. ;  LASA  Travel-Time  Anomalies  for  Various  Epicentral 
Regions;  13  September  1966 
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